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Abstract. Motivated by recent experiments on thin films of the ferromagnetic insulator yttrium-iron garnet 
(YIG), we have developed an efficient microscopic approach to calculate the spin- wave spectra of these 
systems. We model the experimentally relevant magnon band of YIG using an effective quantum Heisenberg 
model on a cubic lattice with ferromagnetic nearest neighbour exchange and long-range dipole-dipole 
interactions. After a bosonization of the spin degrees of freedom via a Holstein-Primakoff transformation 
and a truncation at quadratic order in the bosons, we obtain the spin-wave spectra for experimentally 
relevant parameters without further approximation by numerical diagonalization, using efficient Ewald 
summation techniques to carry out the dipolar sums. We compare our numerical results with two different 
analytic approximations and with predictions based on the phenomenological Landau-Lifshitz equation. 

PACS. 75.fO.Jm Quantized spin models - 75.30.Ds Spin waves - 05.30.Jp Boson systems 



1 Introduction 

In a recent series of experiments 1,2,3,4,5 Demokritov 
and co-workers discovered strong correlations of highly oc- 
cupied magnon states in thin films of the magnetic insula- 
tor yttrium-iron garnet (YIG) with stoichiometric formula 
Y 3 Fe2(Fe04)3. They suggested an interpretation of their 
results in terms of Bose-Einstein condensation of magnons 

, at room temperature. For a proper interpretation of these 
experiments, a peculiar feature of the energy dispersion 
of the relevant magnon band in finite YIG films is impor- 
tant: due to a subtle interplay between finite-size effects, 
short-range exchange interactions, and long-range dipole- 
dipole interactions, Ej- exhibits a local minimum at a finite 
wave- vector fc m i n , for a certain range of orientations of the 
external magnetic field H e relative to the sample. The ex- 
istence of such a dispersion minimum has been predicted 

1 by Kalinikos and Slavin [6j[T] within a phenomenological 
approach based on the Landau-Lifshitz equation. Unfortu- 
nately, such a phenomenological approach does not pro- 
vide a microscopic understanding of correlation effects, 
which might be important to explain some aspects of ex- 
periments probing the non-equilibrium behaviour of the 
magnon gas in YIG P^IBllillS^IWTOini] . This has moti- 
vated us to study this problem within the framework of 
the usual 1 / S-cxpansion for ordered quantum spin sys- 
tems, which is based on the bosonization of an effective 
microscopic Heisenberg model using either the Holstein- 
Primakoff [12] or the Dyson-Maleev transformation [TJl 
[TJ], and the subsequent classification of the interaction 
processes in powers of the small parameter 1/5. 

The 1/S-expansion has been extremely successful to 
understand spin- wave interactions in ordered magnets [HI 
116) . Previously, several authors have used this approach 



to calculate spin-wave spectra in ultrathin ferromagnetic 
films with exchange and dipole-dipole interactions [17,18, 
H91I20], Moreover, interaction effects such as energy shifts 
and damping of spin-waves in thin films have also been 
calculated within the l/S'-expansion [2Tl[2"2~] . However, in 
order to apply this approach to realistic models for ex- 
perimentally relevant YIG films with a thickness of a few 
microns (corresponding to a few thousand lattice spac- 
ings), one has to evaluate numerically rather large dipolar 
sums |23j to set up the secular matrix whose eigenval- 
ues determine the magnon modes, see Eq. (fTH|) below. We 
use here an efficient Ewald summation technique [24] to 
carry out these summations, which enables us to calcu- 
late the spin- wave dispersions of realistic YIG films. Given 
our numerical results, we can assess the validity of various 
analytical approximations such as the uniform mode ap- 
proximation [6, 20, 25J and the lowest eigenmode approxi- 
mation. 

The rest of this paper is organized as follows: In Sec. [2] 
we introduce the effective Heisenberg model which we shall 
use to describe the experimentally relevant magnon band 
in YIG. We set up the l/S'-expansion and derive the secu- 
lar equation which determines the magnon dispersion. In 
Sec. |3] we present our results for the magnon spectra of 
YIG. We first discuss our numerical results, which are ob- 
tained by evaluating the roots of the secular determinant 
without further approximation, using the Ewald summa- 
tion technique described in the appendix to evaluate the 
necessary dipolar sums. We then discuss in Sees. 13.21 and 
13.31 two approximate analytical methods for obtaining the 
dispersion of the lowest magnon band. A comparison with 
our numerical results allows us to estimate the accuracy 
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of these approximations. Finally, in Sec. [4] we present our 
conclusions and give an outlook for further research. 



2 Effective Hamiltonian for YIG films 

Experimental and theoretical research on YIG has a long 
history, as reviewed, for example, in Ref. [26j . Some dis- 
tinct advantages of YIG are that this material can be 
grown in very pure crystals and has a very narrow fer- 
romagnetic resonance line, indicating very low spin-wave 
damping. Actually, YIG is a ferrimagnet at accessible 
magnetic fields and has a rather complicated crystal struc- 
ture with space group Ia3d (see Refs. [2"M2"T] ) and 20 mag- 
netic ions in the primitive cell. Fortunately, on the energy 
scales relevant to experiments [Tll2~, 3 4,5,8,9,1 OlfTT] only 
the lowest magnon band is important, so that we can de- 
scribe the physical properties of YIG at room temperature 
in terms of an effective spin S quantum Heisenberg ferro- 
magnet on a cubic lattice with lattice spacing [27) 



= 12.376 A , 



(1) 



The effective Hamiltonian contains both exchange and 
dipole-dipole interactions, 



H — — — ^ JijSi ■ Sj — [iH e ■ ^ Si 



1 ^ M 2 



E 



2 



3{Si • Rij)(Sj • Rij) Si • Sj 



where the sums are over the sites Ri of the lattice and 



R; 



Rij/\Rij\ 



are unit vectors in the direction of 



Rrj = R4 - Rj = x tj e x + y^By + z y e z . Here, fi = g\i B 
is the magnetic moment associated with the spins, where 
g is the effective (/-factor and [Ib — eh/(2mc) is the Bohr 
magneton. The exchange energies Jy = J(Ri~Rj) decay 
rapidly with distance, so that it is sufficient to include only 
nearest neighbour exchange couplings in Eq. @, setting 
Jij = J if Ri and Rj are nearest neighbours, and Jy = 
otherwise. Note that we neglect surface anisotropies which 
might be present in experiments, especially at the surface 
of the film which is attached to the substrate. Experimen- 
tally, the material YIG is characterized by its saturation 
magnetization [2FJ 



A-kMs = 1750 G, 



(3) 



and the exchange stiffness p C x of long-wavelength spin- 
waves 129 , 



Pcx 



JSa 2 



5.17 x 10" 13 Oem 2 . 



(4) 



If we arbitrarily set the effective g-factor equal to two [55] 
so that /j, — 2/j.b, we obtain from Eqs. (J3J and for the 
effective spin 



s 



14.2, 



li 



(5) 
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Fig. 1. (Color online) Orientation of our coordinate system 
for an infinitely long stripe of width w and thickness d. We 
assume that the external magnetic field H e e z is parallel to the 
long axis (which we call the 2-axis) of the stripe. 



and for the nearest neighbour exchange coupling 



J = 1.29 K. 



(G) 



Our above estimates for S and J differ slightly from the 
values given in Ref. [3D] . Note that the effective spin S « 
14.2 is quite large, so that an expansion in powers of 1/5* is 
justified. Introducing the dipolar tensor D"? = D a ^(Ri — 
Rj), 



3R%R?j - S# 



d 2 



i 



(1 <;/,•;•, R, 



(7) 



we can write our effective Hamiltonian ([2]) in the compact 
form 



H = - 



S?S?-hJ2s*, (8) 



ij afl 



where the z-axis of our coordinate system points into the 
direction defined by the magnetic field H e and we have 
introduced the associated Zeeman energy, 



h = [i\H e 



(9) 



We have thus related the set of parameters a, S, J, h ap- 
pearing in our effective Hamiltonian to experimentally 
measurable quantities. 

To proceed, we restrict ourselves to the description of 
an infinitely long stripe of width w and thickness d = Na, 
consisting of N layers. For the stripe geometry shown in 
Fig. [T] where the magnetic field points in any direction 
parallel to the stripe the classical groundstate is a satu- 
rated ferromagnet. Therefore we can expand the Hamil- 
tonian in terms of bosonic operators describing fluctua- 
tions around the classical groundstate, using either the 
Holstein-Primakoff or Dyson-Maleev transformation [T51 
116] . The resulting bosonized spin Hamiltonian is of the 
form 



H — Hq + H n 



(10) 



It turns out that Holstein-Primakoff and Dyson-Maleev 
transformations give different results for H n with n > 4, 
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but the expressions for H~2 and H3 are identical in both 
transformations. The classical ground state energy is 



J 4 



D zz 
y 



Df 



2S 



and the quadratic part of the Hamiltonian reads 



with 



E 



A l3 b\b 3 



Bi 



bib. 



Aij — Sijh + S(5{j S ' Ji n Jij) 



(11) 



(12) 



with the amplitude factors [21|l22j 
A k (x tj ) = e~ lk ' r A(x t - Xj ,r) , 

r 

= SJ k { Xij ) + Stj [h + Sj2D z Q z (i 

n 



S 



(16a) 



D%*{ Xij ) - 2iD^( Xij ) - I)?;' {■>;,)] , (16b) 



and the exchange matrix 

Jk(xij) = J[$ij{6 - Sji - 5jN 

-2(cos(fcj,a) + cos(k z a))} - 6 ij+1 - <%-i] • (17) 



+ 5 



TiXX 

ij 



D yy 



Bi 



%E D i 



Df; 2;/;;;' t)»;<]. 



(13a) 
(13b) 



Since H n /S 2 = 0{l/S n / 2 ) and the effective S is large, 
we expect accurate results even if we only keep the first 
two terms. In order to proceed we have to keep in mind 
that a stripe with thickness d and width w is obviously 
not translationally invariant, so that we cannot simply 
use a full Fourier transform to diagonalize the Hamilto- 
nian. Because in the experimentally studied samples [TJ[H 
[HlHllMSligifTOirn] the width w of the stripe is much larger 
than the thickness d, we can assume that w is practically 
infinite, so that our system can be considered to have dis- 
crete translational invariance in the y- and z-dircctions. 
We may then partially diagonalize H2 in Eq. (TT2")) by per- 
forming a partial Fourier transformation in the j/z-plane. 
Setting Ri = (a^r,-) with = (yi,Zi), and introducing 
the two-dimensional wave- vector k = (k y ,k z ) in the yz- 
plane, we expand 



h = 



1 



k 



bk{xi) 



(14) 



where N y and N z is the number of lattice sites in y and z 
direction. Our Hamiltonian H2 in Eq. (|12l) takes then the 
form 

H 2 = Y E [ A k( x ij) b U X i) b k(Xj) 



k Xi,Xj 

2 

BUM 



b k {xi)b-k{xj) 
bi{^)bl k { Xl ) 



(15) 



3 Spin-wave spectra of YIG 
3.1 Numerical approach 

If the lattice has N sites in the x-direction, then for 
fixed two-dimensional wave-vector k there are N allowed 
magnon energies E n k, n = 0, . . . , N — 1, which are given 
by the positive zeros of the secular determinant 



det 



( Ek — Afc 



0. 



(18) 



'k — Afc 

Here, the N x iV-matrices and B& are defined by 

[Afc]y = A k (xij), (19a) 



[B fc ] 



B k (xij), 



(19b) 



where 1 < i, j < N label now the lattice sites in the in- 
direction. The condition (jT5J) follows simply from the fact 
that the magnon energies can be identified with the poles 
of the propagators of the Gaussian field theory defined by 
the quadratic Hamiltonian (fTS"]) . Note that for TV = 1 the 
condition (|18[) correctly reduces to the diagonalization of 
the Hamiltonian (|15p via Bogoliubov transformation. To 
obtain the complete magnon spectrum of the thin film 
ferromagnet we have to calculate the dipolar matrices in 
Eq. (fTo]) which leads to the calculation of the following 
dipolar sums for fixed Xij, 



Tij 

— e -i(k v yij+k z z ij ) 



a/3 
ij 



Vi], z a 



5 a ? 



(4 + 4 + 4) 3/2 (4 + 4 + 4) 5/2 



(20) 
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Fig. 2. (Color online) Spin-wave 
dispersion of a YIG film with 
thickness d = 400a « 0.495 ^m 
in a magnetic field of 700 Oe for 
wave-vectors in the plane of the 
film and for different propaga- 
tion angles k = 0°, 30°, 60° 90° 
relative to the magnetic field 
(from top left to bottom right). 
The black solid curves represent 
the results from the numerical 
approach for the lowest modes. 
The other curves are obtained 
from Eq. (|36[) using the approx- 
imation Eq. (|31|l (dashed) and 
Eq. (J4TJ) (dash dotted) for the 
form factor f k . The thick dot- 
ted line labelled Eo marks the 
energy of the ferromagnetic res- 
onance given by Eq. t| 2 X p : the 
thick dotted line labelled E„ is 
the energy E a = h + A/2 of 
the classical surface mode for 
Ok = 90° at large wave-vectors, 
see Eqs. ({22 ] |23 ]> . 



where ^2 excludes the term = Zij = when Xij = 0. 
As these sums are slowly converging and previously used 
summation techniques \22\ are not very efficient for small 
wave- vectors needed for the dipolar-dominated spectrum 
of YIG, we use the Ewald summation technique to obtain 
fast convergence. Details of the calculation can be found 
in the appendix. To determine the dispersion of all modes 
we numerically calculate the roots of Eq. (fTSf which can 
be easily done up to a film thickness d rs 7 /im. Using the 
fact that n/h « 2.803 x lCr 3 GHz/Oe we plot all results 
in units of GHz which is most convenient for experiments 
using microwave resonators and antennas to detect mag- 
netic excitations. Typical magnon dispersions for differ- 
ent angles Ok between the propagation direction and the 
magnetic field are shown in Fig.[2]for a film with thickness 
d = 400a. Since our approach includes all effects of the ex- 
change and the dipolar interaction within linear spin wave 
theory, it reproduces semiclassical approximations based 
on the Landau-Lifshitz equation at long wave-lengths. In 
particular, for k — > the dispersion of the lowest mode 
approaches, independently of the thickness d, the classical 
ferromagnetic resonance energy 

E = ^h{h + A) , (21) 

where we introduced the characteristic magnon energy due 
to dipolar interactions, 

A = Att^M s . (22) 

The energy Eq is indicated as a dotted line in Fig. [51 
The spacing between different magnon modes decreases 
with increasing film thickness d; for example in Fig. [3] 
we show the magnon spectrum of a film with thick- 
ness d = 4040a = 5 ^m, which forms already a quasi- 



continuum for energies close to the ferromagnetic reso- 
nance. Moreover in the regime dk <C 1 one observes that 
the energy of the n-th mode deviates from the first mode 
by an energy AE n = p ex ir 2 n 2 /d 2 which reflects the well 
known quadratic behaviour of ferromagnetic spin waves in 
three dimensions. 

Of particular interest is the lowest magnon mode, 
whose dispersion is illustrated in Fig. |4j For propaga- 
tion directions parallel to the magnetic field (Ok = 0) 
we recover the well known [B] minimum of the lowest 
magnon mode at a finite wave-vector fe ro ; n where large 
magnon densities have been detected by their microwave 
radiance [2[5j. As illustrated in Fig.0 the position of the 
minimum depends on the film thickness as expected and is 
less pronounced for ultrathin films. With increasing prop- 
agation angle Ok the minimum becomes more shallow and 
completely disappears for Ok = 90°. From Figs. and [3] 
it is also obvious that for angles Ok > 45° the lowest 
modes shift upwards and tend to hybridize with higher 
modes. However, in the absence of additional symmetries 
the energy levels never cross. The magnon modes with 
finite group velocities v g (k) = 'VkEk tend upwards and 
form the quasi-continuous surface mode. There are no fur- 
ther hybridizations as soon as the mode energies reach the 
energy 

E s = h + A/2 (23) 
of the classical surface mode. 

3.2 Uniform mode approximation 

The dispersion of the lowest magnon band can be derived 
from an effective in-plane Hamiltonian using various ap- 
proximations. In fact, at long wavelengths the dispersion 
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Fig. 3. (Color online) Spin-wave dispersion of a YIG film 
with thickness d = 4040a = 5fim for wave- vectors parallel to 
the external field H e = 700 Oe for Ok = 0° (top) and for 
Ok = 90° (bottom). 
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Fig. 4. (Color online) Spin-wave dispersion of the lowest mode 
of a YIG film with H e = 1000 Oe and d = 5.1^im obtained from 
the numerical solution of Eq. (|18p . Starting from the minimal 
energy E m i n contour lines with the spacing 20 MHz are shown 
to illustrate the rather flat minimum, whereas for larger ener- 
gies the distance between contour lines is 100 MHz. 



of the lowest magnon band can also be obtained within 
the macroscopic approach based on the Landau-Lifshitz 
equation [6|. In the simplest approximation, we ignore the 
fact that the system is not translationally invariant in the 
x-direction and approximate the corresponding eigenfunc- 
tions by plane waves. The lowest magnon band is then 
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Fig. 5. Momentum fc m i n and energy -E m i n of the minimum 
in the dispersion of the lowest magnon mode for Ok = as 
a function of film thickness for different magnetic fields. From 
top to bottom H e = 2100, 1400, 1050, 700, 350, 175 Oe. 



obtained by replacing the operators b k (xi) in Eq. (fT5|) by 



(24) 



Then Eq. (fl"5j) reduces to 

Hf = j2 U k bib k + ^b k b_ k + ^b{bi k 

k L 



with 



ij 



N ■ 

We parametrize the in-plane wave-vectors as 
k = |fc| (cos6>fc e z + sin6>fc e y ) , 



(25) 

(26a) 
(26b) 

(27) 
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Fig. 6. (Color online) Plot of the form factors Eq. (|3T|) for 
the uniform mode approximation (dashed) and Eq. (|47p for 
the lowest eigenmode approximation (dash dotted). For the 
samples used in the experiments the minimum of the lowest 
mode is at dk > 5 where the lowest eigenmode approximation 
is more accurate. 



carry out the integrations over the two infinite directions 
and obtain for x^j =/= 



D% x ( Xij ) = D k ( Xij ), 
D l V ( x ij) = - sin2 @k D k (xij), 
D k( x ij) = - cos2 °k D k (xij), 

D k V ( x ij ) = sin ®k D k (xi 3 ■ ) f ' ^ 



sign(xij) 



where 



D k ( Xij ) = ^\k\e-M^\ 



(28a) 
(28b) 
(28c) 

(28d) 



(29) 



It should be noted that in the derivation of the dipolar 
tensor (|28aH28d]> we have implicitly assumed that Xij ^ 0. 
This has important consequences when we also replace 
the sum J2x by the integral a^ 1 J^^ 2 dx. To properly 
account for the factor (1 — Sij) in the dipole tensor in 
Eq. ([7]) we therefore exclude a sphere of infinitesimal ra- 



dius around x, 



Vij 



in our integrations, giving 



rise to the dipole matrix elements 



T^XX 

u k — 



Dl V 



Airp 2 



3 



47T/i 



2 r 



sm 2 e k (/„ - 1; 



47r^ 2 



-+cos 2 0fe(/fe-l) 



(30a) 
(30b) 
(30c) 



and D^ y = 0. Here we have introduced the form factor [30] 



i^ = 1 _H? + 0( *V), (31, 



fk 



which is shown as a dashed line in Fig. [6] From Eqs. (|30al 
I30c|) it is immediately obvious that our dipole elements 
satisfy the constraint 



D%*+Df+D%* 



o. 



(32) 



which follows also directly from the definition ([7]) of the 
dipolar tensor. In terms of the above dipole matrix ele- 
ments we can now write the coefficients A k and B k as 



A k = h + JS [4 - 2 cos(kyd) - 2 cos(fc z a)] 



-2^ 



Dl V ) 



Dl V ) 



A 
3"' 



(33a) 
(33b) 



The magnon dispersion E k is now obtained by diagonaliz- 
ing the effective Hamiltonian ((25)) via a Bogoliubov trans- 
formation [15pi6j. resulting in 



Eh = \/Al 



\B k \ 2 = y/(A k -\B k \)(A k + \B k \). (34) 



If we expand the terms involving the exchange interaction 
for small k (which is sufficient for the energy range in 
experiments on thin films) we get 

JS'[4-2cos(fc y a)-2cos(fc z a)] ps JSa 2 k 2 = p cx fe 2 . (35) 

We can then simplify the dispersion to 

E k = 

[h + p^k 2 + A{\ - f k ) sin 2 fe ] [h + p c *k 2 + Af k ], 

(36) 

which is indicated by dashed lines in Figs.[H[3]and[71 Obvi- 
ously, this approximation is in good qualitative agreement 
with the numerical result for the lowest mode as long as 
Ok ^ 45°. Eq. (f3"fj)) has been used in Ref. [3U] to discuss 
the role of magnon-magnon interactions in YIG and has 
recently been rederived by Rezende [5D]. Note, however, 
that in Fig. [3] one clearly sees deviations from the nu- 
merical result at intermediate wave- vectors in the interval 
5 % \k\d < 50. By comparing Fig. [5] with Fig. [3] we con- 
clude that for the samples used in experiments [1,2,3,4,5, 
8.9, iUlfTTj the minimum of the dispersion is exactly in this 
range so that better analytical approximations are needed 
for a more accurate description of the magnon dispersion 
in the vicinity of the minimum. 



3.3 Lowest eigenmode approximation 

To derive the dispersion of the lowest magnon mode more 
systematically, suppose that the ip n k(xi) form (for fixed 
fe) a complete set of orthogonal functions with respect to 
the n-direction, i.e., 

^^nkiZi^rnkiXi) = $nm, (37) 

Xi 

^^nk{XiW nk {x 3 ) = Sij. (38) 

n 

We may then expand the operators b k (xi) in this basis, 
b k {xi) = ^ ^nk(xi)b„k, (39) 



Andreas Kreisel et al.: Microscopic spin-wave theory for yttrium-iron garnet films 



7 



where 



bnk = y^^nh( x i) b k(xj 



(40) 



Oft = 90° 



Let us retain in the expansion (|39p only the n — term, 
b k (xi) w ipok(%i)bok, (41) 



If we choose -0ofc (a^j ) = 1/ V (corresponding to the n = 
term in the plane wave expansion) and identify = &ofc 
we recover Eq. (|24p . However, a truncated expansion in 
plane waves seems not to be a good approximation for a 
system of finite width. To improve on the approximation 
(|24|) it is better to expand in terms of the eigenfunctions 
of the exchange matrix given in Eq. (|17p . 



^ ' Jk(Xjj)lpnk( x j) — ^nk?Pnk( x i)- 



(42) 



For open boundary conditions these are standing waves 
with nodes at x — id/2, i.e., 



1pnk(Xi) 



sin[fc n (a; i + d/2)], 



(43) 



where k n = (n + l)7r/d, n = 0, 1, . . . , N — 1. The approxi- 
mation (|4"T]) then reduces to 



(44) 



&fc(zi) « \l — cos(k Xi)b kl 
with fco = 7r/c?, and 

6 fc = ^2 cos(k Xi)b k {xi). 



(45) 



Substituting Eq. (|4"4")l into Eq. (|15p we obtain again an 
effective Hamiltonian of the form (|25|) , but now with 



A k = — ^cos(k Xi) cos(koXj)A k (xij), 



B k = — ^2 cos(k Xi) cos(k x j)B k {xij). 



(46a) 



(46b) 



Again the summations are replaced by integrations which 
can be carried out analytically and result in the dispersion 
of the form of Eq. (|36[) , but with the form factor now given 

by 



fk = 1 - \kd\ 



|fcd| 3 + |fcd|7T 2 + 27r 2 (l + e-l fcd l) 

(fc 2 d 2 +7T 2 ) 2 



1 - — \k\d + 0(k 2 d 2 ). 



(47) 



In Fig. [S]we compare this form factor with the correspond- 
ing form factor (|3ip obtained within the uniform mode ap- 
proximation. Obviously, there are only small differences: 
the linear coefficient in the Taylor series is different, re- 
sulting in a smaller slope of the dispersion (|36[) at k = 0. 
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Fig. 7. (Color online) Enlarged details of Fig. [3] on a lin- 
ear momentum scale: spin-wave dispersion of a YIG film with 
thickness d — 4040a = 5/im for wave- vectors parallel to the 
external magnetic field H e = 700 Oe (Ok = 0°, top) and for 
wave- vectors perpendicular to the magnetic field (Ok ~ 90°, 
bottom). The solid lines are exact numerical results obtained 
from the solution of Eq. (|18|l . The dashed line is our approxi- 
mate expression (|36|) for the lowest magnon band, using the 
uniform mode approximation (|31[) for the form factor. The 
dashed-dotted line is the lowest magnon band with form factor 
(|47[) given by the lowest eigenmode approximation. 



To estimate the validity of these approximations, we com- 
pare them in Figs. [7] and [5] with our numerically exact 
results of the model ([2]) for experimentally relevant param- 
eters. In Fig. [7] we show the relevant details of Fig. [3] on 
a linear momentum scale; obviously for sufficiently large 
wave- vectors the lowest eigenmode approximation is more 
accurate and practically lies on top of the numerically ex- 
act result. On the other hand, for small wave- vectors the 
uniform mode approximation fits better, as illustrated by 
the lower part of Fig. [7] A more quantitative comparison 
between these two approximations for wave-vectors par- 
allel to the magnetic field and different film thickness is 
shown in Fig. [H] Roughly, for d\k\ < 5 the uniform mode 
approximation fits better, while for for d\k\ > 5 the lowest 
eigenmode approximation gives a better agreement with 
the numerical results. As there are not many states at 
small wave-vectors, we believe that the lowest eigenmode 
approximation is more suitable for a quantitative descrip- 
tion of the magnon dispersion, in particular if one is inter- 
ested in physical effects related to the dispersion minimum 
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Fig. 8. (Color online) Accuracy of our two analytical ap- 
proximations for the lowest magnon mode for a YIG film with 
N = 400 layers (top) and N = 4040 layers (bottom). Dashed 
line: difference between the numerical result and the uniform 
mode approximation; dash-dotted line: difference between the 
numerical result and the lowest eigenmode approximation. 

at fe m in- On the other hand, both analytical approxima- 
tions describe the quasi-continuous surface mode for large 
Ok > 45° rather accurately as long as the wave-vectors 
are sufficiently small so that the dispersion curves tend 
upwards due to the presence of exchange interaction [7] 
and hybridize. 



4 Conclusions and outlook 

In this work we have presented a comprehensive discussion 
of the spin-wave spectra of experimentally relevant sam- 
ples of thin films of the magnetic insulator YIG. Starting 
from an effective spin-5 Heisenberg Hamiltonian with ex- 
change and dipole-dipole interactions on a cubic lattice, 
we have used a truncated Holstein-Primakoff transforma- 
tion to obtain an effective quadratic boson Hamiltonian 
whose eigen-energies can be identified with the magnon 
energies. We have then used numerical methods to calcu- 
late the magnon dispersions for experimentally relevant 
films with thickness corresponding to a few thousand lat- 
tice spacings without further approximations. In order to 
carry out the dipolar sums entering the secular determi- 
nant, the use of efficient Ewald summation techniques was 
necessary. We have also estimated the accuracy of two dif- 
ferent analytical approximations for the lowest magnon 



band: the uniform mode approximation (where the trans- 
verse spatial variation of the eigenmodes is ignored) and 
the lowest eigenmode approximation (where the lowest 
transverse mode is approximated by the lowest eigenmode 
of the exchange matrix). For realistic films, the latter ap- 
proximation is more accurate for wave- vectors in the vicin- 
ity of the dispersion minimum. 

In contrast to the phenomenological approach based 
on the Landau-Lifshitz equation 0(7], in our microscopic 
approach it is straightforward to systematically take into 
account interaction effects. In future work we shall care- 
fully derive the momentum-dependent interaction vertices 
of the bosonized effective in-plane Hamiltonian using the 
lowest eigenmode approximation, which according to our 
investigations in Sec. [3] is more accurate in the experi- 
mentally interesting regime of wave- vectors close to k m j n . 
Although rough estimates of the interactions vertices can 
be found in the literature [30] , more accurate micro- 
scopic calculations which properly take into account the 
momentum-dependence of the vertices are needed in order 
to gain a better understanding of the role of spin- wave in- 
teractions in the experiments p~IT2 ,3.4,5.8.9, 1 0]l 1 1| . For 
example, it would be interesting to know the intrinsic 
damping of the lowest magnon mode for wave- vectors close 
to fc m i n due to magnon-magnon interactions. 
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APPENDIX: DIPOLAR SUMS 



For the calculation of the slowly converging sums in 
Eq. PD|) we introduce the quantity 



e -*(fc»V«+fc«2y) 

(4 + y 2 ,+4) 5/2 ' 



(A.l) 



and get the dipolar sums as derivatives, 



Dl y {x i3 ) 



D%. z {xij ) 



t x vi 



' o 2 


-2*- 




dk 2 v 


' d 2 


-2*- 


dk 2 v 


dk 2 z 


' d 2 


d 2 


dk 2 v 


+ dk 2 z + 




W V h{X 



2xf 



h(xij) , (A.2a) 
I h (xij) , (A.2b) 
I h (xij) , (A.2c) 
(A.2d) 



Note the symmetry = D^ x and the relation D^ v = 

D~f with k — k z e y + k y e z . To evaluate the above expres- 
sions, we consider the cases ^ and iy = separately. 
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A Case ^ 

Using the identity D ^ = £{37=^ ~ \*+9\fM} 

9 

poo Ae^^U? 

J x n -^ 2 e- ax dx = ^- n a- n - 1/2 (2n-l)ll, (A.3) £(r?- ~ 3 4) cos(fe z% > 3 /2(r?e) , 

(A 10a) 

where n > 0, Re a > 0, and introducing the dummy 2 2 

variable e to check the results, we can rewrite d vv (t -\ = — 

^ ffi2 Li 3 ^ e \k+g\ m9) i 



IhM - M 2 - yj- j o dt We -^S^ - 34)008(^^)008(^^)^/2(4.) 



e -(%+*y)e* e -i(fcy2/i;>+fcsZw)_ (A. 4) 



X 



We split the integral in two parts and use Ewald's 
method [UK] 

£e- StH V^ = ^£e-^, (A.5) 
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(A.lOb) 

2 

D l v ( x v) = *-^2- si g(^i) £( fc y + g y )f(p, q) 

g 

4 £ 5/2^2 

+i — -=— Hj ^2 VU sin(fcj,j/y) cos(fc z Zy )^ 3/2 (rf i e) , 
V r 

(A.lOc) 



to transform the lattice sum into a lattice sum m recip- , , , , , , , . , . r- , 

, , 11 ,,■ where we have used the abbreviations q = x^Je and p — 

rocal space, where q — q v e v + q z e z is a reciprocal lattice ,, % 'J v L 

vector This yields v v \k + g|/(2 v / e) and have introduced the function 

f(p, q) = e- 2p9 Erfc(p -q) + e 2p<? Erfc(p + g). (A.ll) 

/ 7rg3 /" ^ 1 / 2 e ~ |fc 4e° |2 ~ x % et For ^ ne simple cubic lattice the components of the recip- 

a 2 7 rocal lattice vectors are c/j, = 27rm, <? 2 = 27m, {m, n} G Z. 



4(«ii) = x^ 2 



9 

5 



where 



, (A.6) 



Note that this result is independent of the variable e which 
shifts the weight from the sums in real space to the sums 
in reciprocal space as e — > and therefore simplifies to the 
sums given in Ref. [25] for e = 0. 



In this case we have to exclude the point r = because 
is the Misra function^ which for v = 3/2 and z = x > there is no self interaction. We therefore introduce the 



dt t v e~ zt (A.7) B Case x {i = 

In this case we ha 

:tionQ which for v — 3/2 and z = x > there is no self in 
can also be written as dummy variable x = ye y + ze z , 

Evaluating the integral in Eq. jX6| for p > and g 2 > Usin S E T <Q we can rewrite Eq. (AJJj) and use the 

transformation on the reciprocal lattice 



using 



dt iV2 e -« *e-P /< = 



e" p ~ 9 ^ a 2 ei 



<Z 2 



9 

We subtract l/|cc| 5 from the sum on the left, which is 



_|- ^_ [e~ 2p9 (l + 2pg) Erfc(p — g) equivalent to removing this first term in the sum over r 

in our dipole sum. In the limit nOwe obtain, 
+ e 2 ™(-l + 2 OT )Erfc(p + <?)], (A.9) 8 e 3 / 2 ^^ a 

4 = ^2^ £ [ e ~ P a-P 2 ) + 2^ 3 Erfc( P )] 
we can calculate the derivatives according to Eqs. (|A.2[) ' g 

and finally get the result for the sums of the dipole inter- 8e 5 / 2 '2 _ _ _ 

actions between spins located in different atomic layers, + ,- [e~ r £ (1 - 2r 2 e) + 2\/7r|r|v / eErfc(|r| y/e)] 

The Misra function is defined in terms of the exponential 8e 5 / 2 
integral Ei„(s) = die-^^-" via v?„(^) = Ei_„(«). — T^jT- ( A - 14 ) 
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After taking the derivatives according Eq. (|A.2|) we note 
that we get for the dipolar sums the limit q — XjiS — > of 
Eq. (|A.10[) and therefore the removing of the origin does 
not make any difference in the calculation of the dipol 
sums except for omitting the term r = in the real space 
sums. 
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